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We study the contribution of v e — > v T — > r — > fi transitions to the wrong-sign muon 
sample of the golden channel of the Neutrino Factory. Muons from tau decays are 
not really a background, since they contain information from the oscillation signal, 
and represent a small fraction of the sample. However, if not properly handled 
they introduce serious systematic error, in particular if the detector/analysis are 
sensitive to muons of low energy. This systematic effect is particularly troublesome 
for large #13 > 1° and prevents the use of the Neutrino Factory as a precision 



facility for large #13 . Such a systematic error disappears if the tau contribution 
to the golden muon sample is taken into account. The fact that the fluxes of the 
Neutrino Factory are exactly calculable permits the knowledge of the tau sample due 
to the v e — > v T oscillation. We then compute the contribution to the muon sample 
arising from this sample in terms of the apparent muon energy. This requires the 
computation of a migration matrix M™ which describes the contributions of the tau 
neutrinos of a given energy Ei, to the muon neutrinos of an apparent energy Ej. 
We demonstrate that applying M^ to the data permits the full correction of the 
otherwise intolerable systematic error. 



1 Introduction 

It has been understood for a long time that muons are the best experimental tool to 
study the yet-unknown parameters of the PMNS matrix [1] (the small angle #13 and 
the Dirac-CP violating phase S) at the Neutrino Factory 

In a Neutrino Factory [2,3], two intense beams of neutrinos are produced by the decay 
of muons which have previously been accelerated to high energy and stored in a ring. 
The so-called golden signature [4] is due to the oscillation v e — > v^, followed by the 
interaction v^ — > \x in a detector. The resulting muon has the opposite sign of the 
lepton circulating in the storage ring, thus the name "wrong-sign" muon [2]. 

The conventional detector considered at the Neutrino Factory is a large iron calorime- 
ter, the so-called Magnetized Iron Detector (MIND) proposed in Ref. [5] and whose 
detailed performance was first studied in Ref. [4]. Those original studies focused on 
achieving a very pure sample of wrong sign muons. To that extent, stringent cuts in 
the muon momentum and transverse momentum with respect to the hadronic jet (Q t ) 
were introduced. As a result, the selection efficiency for neutrinos of energies less than 
10 GeV was very low. 

It was first understood in Ref. [6] that eliminating low-energy neutrino bins from the 
oscillation analysis introduced serious problems in the joint determination of 13 and S, 
namely the appearance of correlations and degeneracies (CAD) [6,7,8,9]. Understanding 
those CAD has been the subject of much theoretical work, both in the framework of 
Neutrino Factories (see, for example, Refs. [10,11,12,13,14,15] and references therein) or 
at other facilities such as Beta-Beams [16,17,18,19,20,21,22] or Super-Beams [23,24,25]. 

Although an iron calorimeter appears limited when it comes to detecting low energy 
muons, refined recent analyses [26] have shown a much better neutrino efficiency than 
the original studies. The ability to include low-energy neutrino bins mitigates greatly 
the CAD problem at the Neutrino Factory [28]. 

The wrong-sign muon sample is contaminated by backgrounds other than the oscillation 
signals. Those backgrounds are reduced by experimental cuts, although a high-efficiency 
analysis for low energy muons must tolerate a larger background, coming primarily 
from neutral currents in which a hadron manages to fake a low energy muon. Those 
backgrounds, however, introduce essentially a "white noise" to the signal (they contain 
no information of the oscillation being studied) and can be statistically subtracted. 

However, the oscillation v e — > v T introduces a small but significant sample of "bona- 
fide" wrong-sign muons, due to the subsequent production of a wrong-sign tau which 
decays in turn in a muon of the same sign. Notice that the wrong sign taus are not 
background but as good signal as the wrong sign muons, but the resulting muons, if 
not properly treated, introduce a systematic error in the muon sample. 

There are two reasons why this systematic error is troublesome: 

(1) In the Neutrino Factory, the total neutrino energy is computed by adding the 



energy of the muon to the energy of the hadronic jet. This operation yields a 
wrong result when the muon comes from a tau decay and it is detected at MIND, 
since no additional information regarding the neutrino missing energy in the decay 
can be provided (in contrast to the case of ECC or Liquid Argon detectors); 
(2) Muons from tau decay tend to accumulate in low-energy muon bins, since the 
missing neutrinos in tau decay result in a "secondary" muon which has, on average 
1/3 of the tau energy. 



On the other hand, the joint determination of #13 and S is particularly sensitive to the 
low-energy muon neutrinos. The small contribution of the taus is in fact significant at 
low energy, in particular for large #13. What is worse is that it is not white noise, since it 
carries oscillation information different from that of the v e — y v^ golden channel. Notice 
that the muons from the tau sample pollute all measurements that use muons in the 
final state as the signal sample. This problem was already discussed in the context 
of precision measurements of the atmospheric parameters (Am|,^3) at a Neutrino 
Factory in Ref. [27]. In that paper, the authors were interested in the study of the v^ 
disappearance transition v^ — y v^ at the Neutrino Factory, with the signal represented 
by the right-sign muon sample in the INO detector. The problem arises, in this case, 
from r's generated through the leading oscillation channel v^ — y is T , with subsequent 
CC interactions v T N . These r's will eventually decay (with a ~ 20% branching ratio) 
into muons, thus contaminating the right-sign muons sample, with a resulting precision 
loss in the (Am^,^) measurement. 



Trying to eliminate muons from wrong- or right-sign taus by using kinematical criteria 
is not a sound procedure, since any cuts which suppress them sufficiently will also 
deplete the low energy muon neutrino bins. Nor is it necessary, since the problem can 
be correctly treated by using migration matrices, a technique already introduced in 
Ref. [17] to correct the effect of computing the neutrino energy from the observed 
lepton in quasielastic events (QE). In this paper, we estimate the systematic error 
introduced in the joint measurement of #13 and 5 if the effect is not treated properly 
and demonstrate that it ruins the sensitivity of the neutrino factory for large #13. 
Next we show that such systematic error disappears if the tau contribution to the 
golden muon sample is taken into account. The fact that the fluxes of the Neutrino 
Factory are exactly calculable permits the knowledge of the tau sample due to the 
v e — y v T oscillation. We then compute the contribution to the muon sample arising 
from this sample in terms of the apparent muon energy. This requires the introduction 
of a migration matrix M^ which describes the contributions of the tau neutrinos of 
a given energy E t , to the muon neutrinos of an apparent energy Ej. We demonstrate 
that applying M^ to the data permits the full correction of the otherwise intolerable 
systematic error. 




w 



10 20 30 40 50 o.o 




Fig. 1. Wrong-sign \x identification efficiency at the Neutrino Factory as a function of the 
reconstructed neutrino energy from Ref. [4] (left) and Ref. [26] (right). 

2 The source of the r-contamination problem 



Figure 1, shows the muon neutrino reconstruction efficiency in the original MIND 
analysis from Ref. [4] (left) and for the improved analysis discussed in Ref. [26] (right). 
Notice that in the early analysis the efficiency below 10 GeV is practically zero. All 
the subsequent studies who used this efficiency curve had, thus, an effective low-energy 
threshold around 10 GeV. 

However, all the studies on the degeneracy problem at the Neutrino Factory have shown 
that the signal content of these low-energy bins is extremely important to solve, or at 
least considerably mitigate the CAD problem: some of the degeneracies are strongly 
energy-dependent, and non- vanishing signal for all ranges of neutrino energies are re- 
quired to distinguish the true solution from its "clones" [13]. For this reason recent 
studies have focused in achieving better neutrino efficiency at the cost of accepting 
higher backgrounds. In Figure 1 (right), the efficiency for events with reconstructed 
neutrino energy in the energy range E u e [5, 10] GeV is around 60% and, even below 
5 GeV, a non-negligible efficiency is found. 

The improved efficiency in the low-energy part of the neutrino spectrum, however, 
has as the drawback that a previously irrelevant background becomes now potentially 
harmful. The problem arises from the contamination of the wrong-sign muon sample 
by wrong-sign muons produced in the decay of wrong-sign r's. It is well known [4,11] 
that v e oscillate into v^ (the "golden" channel) and v T (the "silver" channel) with 
similar rates for the neutrino energies and baselines considered in a standard multi- 
GeV Neutrino Factory setup. Oscillations of v e into v T will produce r's through v T N 
CC interactions within the detector. Approximately 20% of the r's will, eventually, 
decay into muons. Notice that these wrong-sign muons from wrong-sign r's escape 
essentially all filters designed to kill the dominant backgrounds and directly add to the 
wrong-sign muon sample. 
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Fig. 2. The fraction of r contamination of the golden muon sample as a function of the 
reconstructed neutrino energy, for #13 = 2°. Left: L = 4000 Km; Right L = 7500 Km. Top: 
5 = -90°; Bottom: 5 = +90°. 

Certainly if the events were to be divided in muon energy bins (rather than in bins 
of neutrino energy), muons arising from tau decay could be directly treated as an 
additional source of signal. However, the total energy of the neutrino is a fundamental 
input to separate high-energy charged currents from the low-energy dominant neutral 
current background. Thus, in the standard MIND analysis, the neutrino energy is 
induced by adding the energy of the muon and that of the hadronic jet. 

This approach, however, implies that the sample of wrong-sign muons from the decay 
of wrong-sign taus will be distributed erroneously in neutrino energy bins, since the 
missing energy in the r — > v T v il [i~ decay cannot be measured 1 1. The muons from tau 
decay will result, therefore, in a contamination of the wrong-sign muon sample by 
events whose parent neutrino energy is reconstructed wrongly: this is what has been 
called the " r-contamination" problem in [27]. 

In Fig. 2 we show the fraction of muons coming from r-decay that can be found in 
the wrong-sign muon sample after binning in the reconstructed neutrino energy, using 
the new MIND efficiency, for a 25 GeV Neutrino Factory with detectors located at 
L = 4000 Km and at the so-called magic baseline, L = 7500 Km (see Ref. [28] for a 
discussion of why those are the two optimal baselines). The data are shown for #13 = 2° 



1 Notice that this situation is different from what happens in an Emulsion Cloud Chamber, 
where the identification of the r decay vertex allows for a separation of " golden" from " silver" 
muon samples [11]. 



and two values of S, 8 — ±90°. In the left panels, we show the level of r-contamination 
at the L = 4000 Km baseline and in the right, at L — 7500 Km. 

In all cases, a significant r-contamination can be observed below 5 GeV. At the 
shorter baseline — which is optimized to increase sensitivity to 5 — , the amount of r- 
contamination is (^-dependent: we find that 40% (20%) of the muons are produced 
through r-decay for S = —90° (90°). Not so at the long baseline which has been chosen 
to be largely insensitive to S. Here we find that r-contamination below 5 GeV is about 
60%, independently of the value of 5. With the exception of the top left panel (cor- 
responding to # 13 = 2°, 5 = —90°,L = 4000 Km), the r-contamination in the energy 
range E v G [5, 10] GeV is at the percent level. Above 10 GeV, the r-contamination 
becomes negligible (for a 25 GeV Neutrino Factory). For this reason, in previous stud- 
ies of the Neutrino Factory performances obtained using the original MIND efficiency 
(i.e., the muon identification efficiency of Fig. 1, left), the r-contamination problem 
was absent in practice. 

Since the far baseline contributes to data with a substantially <5-blind component (both 
for the true golden channel data and for the corresponding r-contamination), it is not 
able to remove the distortion induced by a wrong treatment of muons from wrong-sign 
taus. For this reason, in the next section we will focus on the intermediate baseline, 
L = 4000 Km, and show only at the end of the section that the combination of the two 
baselines does not solve the problem. 



3 The impact of the r-contamination problem 

In Fig. 3(left) we show the result of a fit to the data in the (#13, 5)-plane performed for 
the input value # 13 = 6.8° (corresponding to the present estimate from the global fit to 
solar, atmospheric and LBL data from Ref. [29]) for three representative values of 5, 
5 = 160°, 30° and —90°. Data have been obtained for a 25 GeV Neutrino Factory, a 50 
Kton MIND located at L = 4000 Km from the source, with 5 x 10 20 useful muon decays 
per year per baseline and 5 years of running time with each muon polarity. Events are 
binned in the reconstructed neutrino energy, with five bins of constant size AE V = 5 
GeV. The input parameters of the simulation, in addition to (#13, 5), have been kept 
fixed to: Am 2 21 = 7.6 x 10~ 5 eV 2 , Am 2 2 = 2.5 x 10~ 3 eV 2 , 9 12 = 33° and # 23 = 42° 
[29]. The data correspond to a MonteCarlo simulation that includes the golden muon 
sample and the r-contamination. On the other hand, the fit to data (in this section) 
is always performed using the theoretical distribution of golden muons (cross-sections 
and efficiencies taken from Fig. 1, right, are properly taken into account). 

It can be clearly seen that the best-fit point does not coincide with the input value 
(represented by the dot), independent of the particular choice of S. At the intermediate 
baseline (left panel), we see that the input value lies generally at the border of the 3<r 
contour of the fit, but for the case of 5 = —90°, when it lies even further away. Notice 
that the x 2 at the best-fit point is, in all cases, rather large: for 5 = 30° (the case in 
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Fig. 3. Left: Ax 2 contours at 1, 2 and 3 a (2 dof's) of a fit of the golden muon theoretical dis- 
tribution to simulated data that include the t- contamination. The data have been produced for 
#13 = 6.8° (corresponding to the present best-fit value from the global fit to solar, atmospheric 
and LBL data, from Ref. [29]) and 5 = 160°, 30° and —90°. The dot represents the input 
value. Right: Test of the hypothesis that a simulation of the data that includes the effect of 
the t- contamination can be fitted with the golden muon theoretical distribution. In the regions 
to the right of the contour lines, the hypothesis can be rejected at 1, 2 or 3a (from left to 
right), assuming the goodness- of -fit statistics follows the x 2 distribution with n = 8 dof's. In 
both panels, data refer to the L = 4000 Km baseline. 

which the distance between the best-fit point and the input value is the smallest), we 
get for the normalized x 2 distribution Xmin ~ 3.63 (for n = 8 dof's). These results 
imply that less than 0.01% of the time we would have got a fit to the data worse than 
what we have found, within the hypothesis that the data are distributed accordingly to 
the golden muons theoretical distribution. In practice, the hypothesis can be rejected 
at more than 3.9a for all of the considered input pairs. 

The result of this analysis can be generalized to different values of the true # 13 . In 
Fig. 3 (right) we present the result of a test of the hypothesis that a simulation of 
the data that includes the effect of the r-contamination can be fitted by a theoretical 
distribution including golden muons, only. The area to the right of the contour lines 
represent the region of the (#13, 5) parameter space for which the hypothesis can be 
rejected at 1, 2 or 3a (from left to right), assuming the goodness-of-fit statistics follows 
the x 2 distribution with n = 8 dof's. The hypothesis that the golden muon theoretical 
distribution can fit the data can be rejected at more than 3cr for a true #13 > 6° for 
almost any value of S (but for relatively small regions around 5 = ±90°). In this part 
of the parameter space, clearly, we would be obliged to use a different distribution to 
fit the data. 

The scenario can, in a sense, be even worse when # 13 is smaller than 5°. In this part of 
the parameter space, we cannot exclude with high statistical significance the hypothesis 
that the data are statistically distributed accordingly to the golden muon theoretical 
distribution, although we would still get a poor fit to the data. Moreover, for input 
parameters in this region the best-fit point of the x 2 fit to the data is always located far 
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Fig. 4. Left: the value of the relative shift in #13 at the best-fit point; Right: the value of the 
absolute shift in 5 at the best-fit point. In all cases, results are shown as a function of the true 
5. Different lines correspond to different true #13 in the range #13 £ [3°, 10°]. Data have been 
obtained at L = 4000 Km. 

from the input value. The distance between the best-fit point of the x 2 fit to the data 
obtained using the golden muon theoretical distribution and the input values (6*13,5) 
is shown in Fig. 4 for several values of the true #13 in the range #13 G [1°,5°], as a 
function of the true 5. Since the considered range for sin 2 2#i3 runs over several order 
of magnitude, in Fig. 4 (left) we present the relative error in #13 (i.e., the distance in 
#13 between the best-fit point of the \ 2 distribution and the input value normalized to 
the input value, A813/813). On the other hand, since at present S can assume any value 
in the range [0, 2ir[, the distance in S is shown in Fig. 4(right) in terms of the absolute 
error, AS. 

As can be seen from Fig. 4(left), the relative error on 813 can be as large as 17%. The 
maximal relative error is found in the region of negative S and small #13, #13 G [1°, 2°] 
(sin 2 2#i3 G [1,4] x 10~ 3 ). The absolute error in the measurement of the CP- violating 
phase can also be huge: from Fig. 4 (right) we can see that \AS\ can be as large as 
40°. The maximal error is found, again, in the region of negative S, but for relatively 
large 813. For the same regions of the parameter space we have found that, combining 
the intermediate and the magic baseline data, the #13 relative error and the S absolute 
error can still be as large as 4% and 12°, respectively. Clearly, such large errors prevent 
the high-energy Neutrino Factory setup from measuring the value of the mixing angle 
#13 with precision at the level of the percent. If #13 is smaller than #13 ~ 5°, a wrong 
treatment of the r-contamination will prevent the use of the Neutrino Factory as a 
precision machine and will introduce a systematic error in the joint measurement of 
6*i3 and S. 

It is worth noting that the r contamination problem could also affect other observables, 
like the CP-discovery potential, that we define as the capability of a Neutrino Factory to 
measure a value of the CP phase 5 different from the CP-conserving cases S = 0, ±71", at 
some confidence level 2 1. A loss in sensitivity could arise when the silver muons, properly 
taken into account, add to the golden sample: in that case, due to the fact that the 



We thank Thomas Schwetz for rising this point. 



CP violating terms in the v e — > z/ M and v e —¥ v T probabilities have opposite signs, we 
would expect a suppression of the bulk of the events which depend on the CP phase 
S. We carefully checked that this is not the case, mainly because the silver statistics is 
irrelevant at very small #13 when compared with the golden channel. 

In this section we have shown at length that a wrong treatment of the r-contamination 
is extremely troublesome if the value of #13 is larger than 1° (as suggested by present 
three-family oscillation fits) and prevents the use of the Neutrino Factory as a precision 
facility for large 13 . Similar results on the impact of r-contamination in the atmospheric 
sector have been found in Ref. [27], where it was shown that the measurement of 
#23 when r-contamination is taken into account has an error that is almost twice as 
large as when muons from r's are correctly removed at 99% CL for the input pair 
#23 = 42°, Am 2 ^ = 2.4 x 10~ 3 eV 2 . This means that the precision on 623 falls from 5% 
to ~10%. As a consequence, the r-contamination severely reduces the ability of the 
high-energy Neutrino Factory setup to discriminate a non-maximal 823 from #23 = 45° 
and the capability of the Neutrino Factory to solve the octant degeneracy [30]. On the 
other hand, the measurement of the atmospheric mass difference Am^ seems to be 
less affected, and it remains at the percent level. 

Notice that one of the main difference between Ref. [27] and our analysis, apart from 
the study of different channels and parameters, is that in Ref. [27] the signal is studied 
as a function of the final muon energy, whereas in this paper we analyse data as a 
function of the reconstructed neutrino energy. Although binning in the final muon 
energy allows the addition of the r-contamination as a signal source to the wrong-sign 
muons sample, as we have explained in Sect. 2 this prevents an optimal use of the 
detector hadronic calorimetry that, in turn, translates into a worse treatment of the 
backgrounds. A detailed analysis of the r-contamination impact on the measurement 
of the atmospheric parameters 623 and Am^ when data are analysed as a function of 
the reconstructed neutrino energy can be easily done, but it is beyond the scope of the 
present paper and it will be presented elsewhere. 



4 The solution to the r-contamination problem 

In the previous sections we have discussed at length the problems related to fitting, 
using the theoretical golden muon distribution, simulated data that include the true 
golden muon sample and its corresponding r-contamination. We will now explain how 
the simulation is performed and how we can use the same procedure used to produce 
the data to actually solve the problem. 

Consider a v T of energy E„ T , interacting in MIND and producing a wrong-sign r of 
energy E T together with a hadronic jet of energy E^. The heavy lepton decays subse- 
quently (~ 20% of the times) into two neutrinos and a muon, with E T = E^ + E miss , 
where E miss is the missing energy carried away by the two neutrinos in the r-decay. 



We have, therefore, 

E Ut = E T + E h = (E^ + E rniss ) + E h . (1) 

Experimentally, we observe the secondary muon and a hadronic jet, a signal essentially 
indistinguishable from that of a wrong-sign muon from CC v^ interactions. However, 
in this latter case, the addition of the (primary) muon energy E^ and of the hadronic 
jet energy Eh results in the correct parent v^ energy, E v = E^ + Eh- On the other 
hand, in the former case the addition of the (secondary) muon energy E^ and of the 
hadronic jet energy Eh results in the wrongly reconstructed fake neutrino energy E fake- 
As can be immediately seen from eq. (1), the relation between E Vt and Ef a j, e is simply: 

Efake = E^ + E h = E Vt — E miss . (2) 

The decays of the produced r would result in muons with all the possible energies 
between m M (the case in which the neutrinos carry most of the energy of the r) and the 
r energy (the case in which the neutrinos carry no energy). If we divide the continuous 
distribution of the r three-body decay in discrete fake neutrino energy bins, we find that 
for a monochromatic v T beam of energy E Ut , the final muon will be assigned to a given 
fake neutrino energy bin of energy Ej with probability Vj(E Vt ), where j = 1, . . . , N£ n . 
The neutrino factory does not produce a monochromatic v T beam, however. Since the 
beam's flux is exactly calculable, we can compute the distribution of v T of a given 
energy E Vt and divide them into v r energy bins of energy EJ, where % = 1, . . . , N^ in . 
The ensemble of the probability vectors Vj(EJ), for % and j running over all the v^ and 
v T energy bins, is represented by the migration matrix M^. 

Mij is, of course, not measured by experimental data taking. However, knowing the NF 
neutrino flux, the oscillation probabilities, the differential v T N cross-section and the 
differential r — > fx decay width, we can compute statistically My through a MonteCarlo 
simulation of the events. This is the procedure we have used to simulate data that 
include both the true golden muon signal and the corresponding r-contamination. The 
M^ migration matrix has been computed using the GENIE neutrino generator [31] 
with 10 6 simulated z/ r 's per neutrino energy bin and 25 bins in the range E Vr e [0, 25] 
GeV. The explicit form of the migration matrix M^ is depicted in Fig. 5. In this figure, 
we show the statistical distribution of the fraction of the events with v T of energy E\ 
that produce a wrong-sign muon whose energy, combined with the hadronic energy Eh, 
will be erroneously assigned to the v^ energy bin Ej. 

After having computed M^ , the number of wrong-sign muons in a given neutrino energy 
bin is 



Ni= E 

i=l,N biv 



j=l,N bin 



(3) 

Once we know the theoretical distribution of the expected experimental muon sample, 

10 




Fig. 5. The migration matrix Mij. 

including both the true golden muon component and the corresponding r-contamination, 
we can use eq. (3) to fit the experimental data. The results of this fit for the same input 
parameters as in Fig. 3 (left) (i.e., #13 = 6.8° and 5 = 160°, 30° and —90° from top to 
bottom) are shown in Fig. 6. As before, simulated data have been obtained for a 25 
GeV Neutrino Factory, a 50 Kton MIND located at L — 4000 Km from the source, 
with 5 x 10 20 useful muon decay per year per baseline and 5 years of running time 
with each muon polarity. Events are binned in the reconstructed neutrino energy, with 
five bins of constant size AE U = 5 GeV. The other parameters used in the simula- 
tion are: Am 2 21 = 7.6 x 10" 5 eV 2 , Am 2 2 = 2.5 x 10~ 3 eV 2 , 12 = 33° and 6 23 = 42° 
[29]. As it can be seen, the r-contamination problem has been completely solved and 
the best fit points coincide with the input points for all considered input pairs. The 
dashed lines have been obtained considering the golden muons only whereas the solid 
ones take into account the contribution of the silver events. A part from a moderate 
5 dependence, which does not produce any relevant precision modification, we cannot 
appreciate any major differences among the two sets of data. We also checked that this 
situation is common to other points in the (#i3, <5)-plane and we safely conclude that 
the silver sample does not have any significant impact on the determination on the two 
unknowns (as also found for the CP discovery potential). 



5 Conclusions 



In this paper we have shown that: 
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Fig. 6. Ax 2 contours at 1, 2 and 3 a (2 dof's) of a fit of the corrected muon theoretical 
distribution to simulated data. Dashed lines have been obtained considering the golden muons 
only whereas the solid ones take into account the contribution of the silver events. The data 
have been produced for #13 = 6.8° (corresponding to the present best-fit value from the global 
fit to solar, atmospheric and LBL data, from Ref. [29]) and 5 = 160°, 30° and —90°. The dot 
represents the input value. Data refer to the L = 4000 Km baseline. 

• The "tau contamination" problem, ignored in early analyses of the Neutrino Fac- 
tory sensitivity introduces, if not properly treated an intolerable systematic error, 
in particular for large #13. In particular, for #13 > 5°, the golden muon theoretical 
distribution is not able to fit the r-contaminated data. Data, in this case, can only be 
properly studied using the final muon energy distribution, thus reducing the capabil- 
ity of the MIND detector to handle the background and keeping it at an acceptable 
level. For #13 G [1°,5°], the error in the joint measurement of #13 and S can be so 
large that it could actually prevent the use of the Neutrino Factory as a precision 
facility. 

• This problem arises because recent analyses of the Neutrino Factory use muon neu- 
trino energies above 1 GeV, while early analysis had an effective threshold of 10 
GeV. It was shown that, as a consequence of this, the CAD problem at high-energy 
Neutrino Factory setups was extremely severe. The inclusion of neutrino energy bins 
below 10 GeV is very important to mitigate the CAD problem, and for this reason 
new (softer) kinematical cuts were applied to the golden muon sample in order to 
recover efficiency at low energy. 

• Taus contaminate low energy bins and cannot be separated from the wrong sign muon 
signal by means of kinematical cuts. However, the use of a migration matrix allows 
to compute their contribution to the signal bin-by-bin. When such contribution is 
properly introduced in the fit, the large systematic error in the determination of 
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#13 and 5 introduced otherwise is solved. Calculation of the migration matrix is 
straightforward since the relative weights of the different v flavours are exactly known 
at the Neutrino Factory. 
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